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We present additional studies of the kinematics of the anomalous events 
observed by the CDF experiment. These events contain a high-py lepton, 
large transverse missing energy (J5 T ), and 2 or 3 high-ii^ jets, one of which 
contains both a displaced secondary vertex and a soft lepton. Previous articles 
detailed the selection and the kinematic properties of these events. In the 
present paper, we use several phenomenological approaches to model these 
data and to estimate their production cross section. 
PACS number(s): 13.85.Qk, 13.20.He, 14.80.Ly 



I. INTRODUCTION 

The CDF experiment has observed |1] an excess of events in the W + 2 and W + 3 
jet topologies in which the presumed heavy-flavor jet contains a lepton in addition to a 
secondary vertex. The rate of these events (13 observed) is larger than what is predicted by 
the simulation of known standard model (SM) processes (4.4 ± 0.6 events expected, including 
single and pair production of top quarks). The kinematical properties of these events are not 
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consistent with what is expected if the excess were due to a statistical fluctuation of the SM 
contributions . This small sample of events has been isolated in a detailed examination 
of the heavy-flavor content of the W+ jet data sample collected by the CDF experiment 
during the 1992-1995 collider run at the Fermilab Tevatron. Jets with heavy-flavor hadrons 
are identified using the SECVTX algorithm that reconstructs secondary vertices (SECVTX 
tags), the jet-probability (JPB) algorithm |§, and the SLT algorithm that identifies leptons 
produced in the decay of b and c-hadrons (SLT tags). Jets containing both a SECVTX 
and SLT tag (supertag) are referred to as superjets. The data set, the analysis tools and 
the Monte Carlo simulations used in this study are the same as those described in Sec. IV 
of Ref. IJ. The characteristics of the 13 events with a superjet are listed in Tables XVI 
and XVII of Ref. 0. 

We are not aware of any model for new physics which incorporates the production and 
decay properties necessary to explain all features of the 13 events with a superjet. In the 
absence of a suitable theoretical model, in this paper we investigate several hypotheses 
with the purpose of evaluating the production cross section of these anomalous events. For 
this purpose, we use tools, names and ideas developed for specific models, none of which 
can be supported or excluded at the current time because of the limited statistics of the 
data. In Sec. [IJ we investigate the hypothesis that superjets are due to the production and 
weak decay of a light spin-0 quark (scalar quark, q). In order to estimate the size of the 
anomalous contribution we fit the yields of W+ jet events with different types of tags with 
a SM simulation implemented with the pair production of a scalar quark and a b quark 
with large invariant mass (~ 220 GeV/c 2 ). The fit result shows that the observed tagging 
rates can be reasonably described with this simple addition, and provides a first indication 
of how much the CDF measurement [|J of the ti production cross section is impacted by the 
presence of events with a superjet. In Sec. |ITJ we use several simulations based upon the 
salient features of the data in order to evaluate the detector acceptance for events with a 
superjet and estimate their production cross section. Section [TV] summarizes our conclusions. 
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II. ESTIMATE OF THE SIZE OF THE ANOMALOUS PRODUCTION 



An interpretation in terms of new physics of the superjet properties, which are described 
in Sec. IX of Ref. Jl], would require the production of a low-mass, strongly interacting object, 
decaying semileptonically with a branching ratio close to 1 and with a lifetime of the order 
of a picosecond. The low mass is required by the fact that the SECVTX and SLT tags are 
contained in a cone of radius 0.4 in the 7] — <j) space around the jet axis. The request that 
the object is strongly interacting is motivated by the fact that several tracks in the superjet 
point directly to the primary vertex, suggesting an emission of gluons during the evolution 
prior to the decay. The large semileptonic branching ratio, finally, is required to explain the 
large fraction of soft lepton tags in jets tagged by SECVTX. 

Since there are no experimental limits on the existence of a charge-1/3 scalar quark 
with mass smaller than 7.4 GeV/c 2 the supersymmetric partner of the bottom quark 

is a potential candidate. In order to explain the observed rate of supertags, one possible 
assumption is that this scalar quark has three-body decay modes to quarks mediated by 
heavy gauge fermions x { m x — m l) anc ^ that X ~ * & i s the only allowed decay. As an 
example, in this study we use the ansatz that superjets are due to the production and decay 
into lev of a scalar quark b with a mass of 3.6 GeV/c 2 and a lifetime of 1.0 ps; the scalar 
neutrino is assumed to be massless. Such a scalar quark is not ruled out by the CLEO 
collaboration || if its decay matrix element results in a lepton spectrum softer than three- 
body phase space, as is the case for the matrix element used by us (see Appendix A). Even 
if such a b quark were ruled out, we would retain this simple model, which fits the data 
reasonably well, as a working hypothesis to estimate the production cross section of the 
anomalous events. In fact, many slightly different conjectures (m? in the range of a few 
GeV/c 2 and 2 < m~ b — m c — < 5 GeV/c 2 ) are still not excluded by any data analysis 
and, because of the limited statistics of the CDF data and the large transverse energy of the 
superjets (~ 60 GeV), they also provide a reasonable modeling of the superjet properties |§. 
We also note that we assume the presence of a c quark in the b decay but this assumption 
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is not required to fit any of the superjet properties. 

The estimate of the production cross section of the events with a superjet requires the 
knowledge of the tagging efficiencies and the detector acceptance, which in turn depend on 
the kinematics of the data. For example, the acceptance is sensitive to the pseudo-rapidity 
and transverse momentum distributions, which are anomalous in the data. We use three 
techniques of increasing complexity to estimate the sensitivity to the details of modeling the 
kinematics. In the first method, we ignore entirely the kinematics of the primary lepton and 
the missing energy, and we model the jet kinematics with the fictitious production of a heavy 
state N of mass 220 GeV/c 2 which decays into bb according to phase space. This mechanism 
is chosen to reproduce the hard E T spectrum of the superjet and the other jets in the final 
state, as well as their invariant mass distribution. The observed rates of tagged events are 
fitted with a SM simulation implemented with the additional production pp — * iV — > bb 
(described in Appendix A). Since this model is the simplest, we describe the fit in some 
detail in subsection A. Using the result of this fit, subsections B-D show that the addition 
of this hypothetical process improves the simulation of other peculiar features of the events 
with a superjet. 

A. Fit of the tagging rates 

We fit the rates of W + n jet events (n=l,4) before and after tagging with the SM 
simulation implemented with the N production using a likelihood technique. In order to use 
all the information contained in the data, we fit nine (k = 1, 9) different classes of tagged 
events: events with one and two SECVTX tags; events with one and two JPB tags; events 
with one and two SLT tags; events with one supertag; events with one supertag and an 
additional SECVTX tag; and events with one SLT and one SECVTX tag in two different 
jets. In total we fit 31 data points. 

The ti cross section and the number of events due to N production are unconstrained 
parameters of the fit, which has 29 degrees of freedom. The remaining eleven contributions 
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Xj to the tagging rates [SECVTX, JPB and SLT mistags, WW, WZ, ZZ, single top, Wc, 
Wcc, Wbb and unidentified- Z with heavy flavors] are also fit parameters. However, we 
constrain each contribution to its expected value, Xj, and within its estimated uncertainty, 
a Xi , using in the log-likelihood function the term 

{Xi - xf 



G , 



2 4, 



The detector acceptance, the luminosity and the values of the SECVTX, JPB and SLT 
tagging efficiencies are five additional fit parameters, which are also constrained to their 
measured value within their uncertainties. 

In summary the fit minimizes the function 

4 9 16 

-ln£ = - ^\nP nk (N nk \p, nk ) -^Gi 

n=l k=l i=l 

where P nk is the Poisson probability of observing N nk events when ft nk events are expected. 

The minimum value of —InC returned by the fit is -1854.7. If each data point were 
equal to the best fit result, we would have found —lnC = -1868.7. The difference be- 
tween —InC and —ln£ corresponds to a x 2 increase of 28 units for 29 degrees of freedom. 
Rates of observed and fitted W+ > 1 jet events before and after tagging are compared in 



Tables| to|m. 

The value of the it cross section returned by the fit is 4.0 ± 1.5 pb, in agreement with 
theoretical expectations |T(J and the D0 result [|TTJ. This cross section corresponds to the 
presence of 41.9 ± 15.7 top events in the W+ > 1 jet sample before tagging. In the same 
W+ > 1 jet sample, the fit returns 52.8 ± 22.1 events due to the production of a massive 
bb pair. We quote numbers of events and not a cross section since the lepton and missing 
energy kinematics have been ignored. 
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TABLE I. Composition of the W+ > 1 jet sample before tagging. This composition is deter- 
mined by a fit to the tagging rates of the data (see text). 



Source 


JU + ljet 


W + 2]et 


JU±3jet 


W+ > 4 jet 


Data 


9454 


1370 


198 


54 


Non-W 


560.0 ±15.1 


71.2 ± 2.6 


12.3 ± 2.0 


5.1 ± 1.6 


WW, wz, zz 


35.8 ± 7.6 


36.2 ± 7.7 


6.2 ± 1.3 


0.9± 0.2 


Unidentified- Z + jets 


234.8 ± 14.5 


38.5± 5.9 


7.9 ± 2.4 


0.7 ± 0.7 


Single top 


14.0 ± 2.2 


7.8 ± 1.3 


1.7 ± 0.3 


0.3± 0.1 


ti 


1.4 ± 0.5 


7.9 ± 3.0 


15.9 ± 6.0 


16.7± 6.3 


W + jets without h.f. 


7955.8± 73.6 


1004.5± 20.5 


110.6 ±6.6 


19.1± 4.9 


Wc 


399.2 ±89.1 


82.1 ±18.3 


10.0 ± 2.2 


1.8±0.4 


Wcc 


171.6 ±36.9 


60.1 ±12.9 


10.4 ± 2.2 


2.2±0.5 


Wbb 


67.6 ± 7.8 


28.8 ± 3.3 


5.1 ± 0.6 


1.4 ±0.2 


N 


2.2 ± 0.9 


28.2 ±11.8 


16.7 ± 7.0 


5.6 ± 2.3 
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TABLE II. Summary of observed and fitted number of W events with one (ST) and two (DT) 
tags of the same kind (SECVTX, JPB, or SLT). "Other" is the sum of all processes other than ti 
and the fictitious N production (see text). 



Source 



Other (ST) 
Other (DT) 
it (ST) 
tt (DT) 
N (ST) 
N (DT) 



W + ljet 



63.68 ±3.89 



0.42 ± 0.12 



0.21 ± 0.06 



W + 2jct 



SECVTX tags 
25.60 ± 1.64 
1.56 ±0.23 
2.59 ± 0.74 
0.58 ± 0.17 
9.32 ± 2.75 
3.08 ± 0.91 



W + 3}et 



5.73 ± 0.39 
0.30 ± 0.04 
5.25 ± 1.50 
2.22 ± 0.63 
6.20 ± 1.83 
1.71 ± 0.51 



W+ > 4 jet 



1.45 ± 0.13 
0.07 ± 0.01 
5.75 ± 1.65 
3.05 ± 0.87 
2.01 ± 0.59 
0.70 ± 0.21 



Fit (ST) 
Fit (DT) 
Data (ST) 
Data (DT) 



Other (ST) 
Other (DT) 
tt (ST) 
tt (DT) 
N (ST) 
N (DT) 



64.30 ± 3.90 
0.00 ± 0.00 
66 




128.46 ± 9.36 



0.41 ± 0.12 



0.24 ± 0.08 



37.50 ± 2.41 
5.22 ± 0.81 
35 
5 



JPB tags 
45.77 ± 3.46 
1.54 ± 0.23 
2.46 ± 0.73 
0.56 ± 0.17 
9.60 ± 2.97 
3.20 ± 1.00 



17.17 ± 1.50 
4.23 ± 0.54 
10 
6 



9.61 ± 0.78 
0.28 ±0.04 
5.12 ± 1.51 
2.02 ± 0.60 
5.70 ± 1.76 
1.93 ± 0.60 



9.21 ± 1.40 
3.81 ± 0.77 

11 

2 



2.06 ± 0.16 
0.06 ±0.01 
5.46 ± 1.61 
2.81 ± 0.83 
1.79 ± 0.55 
0.75 ± 0.23 



Fit (ST) 
Fit (DT) 
Data (ST) 
Data (DT) 



Other (ST) 
Other (DT) 
tt (ST) 
tt (DT) 
N (ST) 
N (DT) 



129.12 ± 9.44 



125 



140.55 ± 7.47 



0.10 ± 0.03 



0.10 ± 0.03 



57.83 ± 5.20 
5.30 ± 0.94 
62 
6 



SLT tags 
45.41 ± 2.38 
0.08 ± 0.01 
1.02 ± 0.32 
0.03 ± 0.01 
8.07 ± 2.43 
0.74 ± 0.22 



20.43 ± 2.66 
4.23 ± 0.63 
21 



10.13 ± 0.55 
0.01 ± 0.00 
2.16 ± 0.68 
0.10 ± 0.03 
5.20 ± 1.56 
0.41 ± 0.12 



9.31 ± 1.82 
3.62 ± 0.76 

12 
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3.72 ± 0.22 
0.00 ± 0.00 
2.55 ± 0.80 
0.13± 0.04 
1.75 ± 0.53 
0.19 ± 0.06 



Fit (ST) 
Fit (DT) 
Data (ST) 
Data (DT) 



140.75 ± 7.47 



146 



54.51 ± 2.99 
0.86 ± 0.32 
56 




17.49 ± 1.44 
0.52 ± 0.17 
17 




8.02 ± 0.77 
0.32 ± 0.07 
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TABLE III. Summary of observed and fitted number of W events with one supertag (ST) and 
one supertag with an additional SECVTX tag (DT) (top table) and with one SLT and one SECVTX 
tag in two different jets (bottom table). 



bource 


TT/ | 1 '4- 

W + ijet 


TT/" i O 

kk + zjet 


W + ojet 


W + > 4jet 


Other (ST) 


3.77 ± 0.35 


1.93 ± 0.18 


0.54 ± 0.05 


0.05 it 0.00 


Utner (L) 1 ) 




n i ^ _i_ n no 
U.lb ± U.Uz 


n HQ _i_ n nn 
U.Uo ± U.UU 


n nn _i_ n nn 
U.UU it U.UU 


tt l ) 


n nQ 4- n nn 
U.Uo ± U.UU 


n on _i_ n iiq 
u.zy ± U.Uo 


n ^7 4- n n^ 
U.0< ± U.Uo 


n 7n 4- n n7 
U. / U ± U.U ( 


tt yu 1 1 




U.U f It U.U1 


n O/i x n no 

U.Z4 It U.UZ 


n Q7 4- n HQ 
U.Of it U.Uo 


AT (Q r P\ 

N (SI) 


n no _i_ n nn 
U.Uz ± U.UU 


1.04 ± U.14 


1 no i n i n 

l.Uz ± U.IU 


n oq 4_ n hq 
U.zo it U.Uo 






i no _i_ n i n 
i.uo in u.iu 


u.uo in u.uu 


n oq _|_ n nQ 

U.ZO It U.UO 


Fit (ST) 


3 82 ± f) 36 


3 76 ± 35 

>.) . t \J _1_ V . >.} 'J 


2 13 it 20 


1 03 it 10 

1 .\ )> ) _1_ \J . A.\J 


Fit (DT) 




1 31 + 12 


90 it 08 


65 it 06 


Data (^T\ 


i 


u 


o 


O 

At 


Data (DT) 




2 


3 





Source 


W + ljet 


W + 2jet 


W + 3jet 


> 4 jet 


Other 




1.10 ± 0.10 


0.44 ± 0.04 


0.20 it 0.02 


ti 




0.22 ± 0.02 


0.60 it 0.04 


0.92± 0.06 


N 




1.17 ± 0.11 


0.89 it 0.08 


0.32± 0.03 


Fit 




2.49 ± 0.23 


1.93 ± 0.18 


1.44 ± 0.14 


Data 




1 





1 
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B. Choice of the N mass 



The additional production iV — > bb is a simple model capable of reproducing the observed 
yield of tagged events as a function of the jet multiplicity (see Tables [TT] and PTI| ). Since the 
tagging rates are not very sensitive to the N mass, we have chosen its value after comparing 
the invariant mass distribution of jet pairs in the data and in the simulation. For this 
purpose, we have added to the original 13 events the four additional events with a superjet 
listed in Tables XVIII and XIX of Ref. ]![]. Of these 17 events, 11 events contain only two 
jets. In the remaining 6 events with three jets, when possible, we combine the superjet with 
the other jet tagged by SECVTX (3 events). If no additional jets are tagged, we select the 
one with the highest transverse energy (3 events). In simulated events due to N production 
containing two untagged jets in addition to the superjet, the highest-^ jet is produced by 
a parton from the iV-decay in 77% of the cases. 

Using the normalization provided by the fit in Section [II A| , the simulation predicts 
10.2 events with a supertag: 6.0 ± 0.6 are due to N production, 1.3 ±0.1 are due to 
ti production and 2.9 ± 0.3 are due to the remaining processes. Figure [l| compares the 
observed and predicted invariant mass distributions. A Kolmogorov-Smirnov (K-S) test of 
these distributions, described in detail in Section VII C of Ref. yields a distance 5° = 
0.37 and a probability P = 10.3%. The result of the K-S comparison of the di-jet invariant 
mass distributions in the data and in the simulation for several Mn values is also shown 
in Figure [l] (the probability P of the K-S test has been converted into a x 2 P er degree 
of freedom). In conclusion, the kinematics of the data is compatible with the additional 
production of a pair of bb partons with an average invariant mass of about 220 GeV/c 2 . 
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50 100 150 200 250 300 

M N (GeV/c 2 ) 

FIG. 1. Di-jet invariant mass distribution in the 17 events with a superjet (•) and in the simu- 
lation implemented with the production of a state N — > bb with a 220 GeV/c 2 mass (histogram); 
the SM contribution is represented by the shaded histogram. The inset shows the yield of the 
reduced x 2 °f the comparison between the di-jet invariant mass distributions in the data and in 
the simulation as function of Mjy. 

C. Transverse momentum of the soft lepton tags 

As discussed in Section IX B of Ref . , the distribution of the transverse momentum of 
the soft lepton tag, p^ LT , in the original 13 events with a superjet is quite anomalous, and 
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has a P = 0.09% probability of being consistent with the SM prediction. Figure |2](a) shows 
the Pt LT distribution in 15 events with a superjet (we add the two events recovered using 
primary plug electrons listed in Table XIX of Ref. [|IJ but not the two events triggered by 
the soft muons because of the trigger py-cut). The probability that the p^ LT distribution for 
these 15 events is consistent with the SM simulation remains small (P = 0.1%). The K-S 
comparison of the p§, LT distribution to a SM simulation implemented with the N production 
(Figure @(a)) yields a distance 5° = 0.34 and a probability P = 11.8%. The improved 
agreement obtained by adding the iV-production can be understood from the comparison 
of the pj> LT distributions in SM processes and in events due to the N production in which 
most of the soft leptons are produced by b decays (Figure @(b)). 




30 40 
p T (GeV/c) 



20 
SLT 



50 



30 40 
p T - (GeV/c) 

FIG. 2. Distributions (a) of the transverse momentum of soft lepton tags in the 15 events with 
a superjet and in simulated events which include the N — > bb production. The SM and the N 
contributions, normalized to the same number of events, are also shown (b) separately. 
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D. Charge correlation between the primary lepton and the soft lepton in the 

superjets 



In 13 out of 17 events with a superjet, listed in Tables XVI to XIX of Ref. [JJ, the charges 
of the primary lepton and the soft lepton tag are opposite. SM processes which contribute to 
events with a superjet do not produce a visible charge correlation; the binomial probability 
of an equal or larger statistical fluctuation is 2.4%. 

In the N —> bb simulation most superjets are produced by semileptonic 6-decays. In 67% 
of the cases, the soft lepton charge has the same sign of the b charge. According to the fit 



in Section [II Aj , 4.2 of the 17 events are attributed to SM processes which do not produce 
a visible charge correlation. If the remaining 12.8 events with a superjet represent the final 
state (bb) + l + + K T (and its charge-conjugate), in which the charges of the primary lepton 
and the b quark are opposite, then in 17 events we expect to find 10.7 events in which the 
primary and soft lepton have opposite charge. The probability that the data are generated 



according to this hypothesis increases from 2.4% to 18. ( 



/o. 



III. ESTIMATE OF THE ACCEPTANCE FOR EVENTS WITH A SUPERJET 



The fit in the Section [II Aj returns 52.8 ± 22.1 events not attributed to SM processes. 



In order to turn this number into a cross section we try two phenomenological approaches 
to estimate the detector acceptance. We first use a W+ Higgs boson simulation where we 
identify the Higgs with the state N which decays to bb. As detailed in Appendix B, we 
arbitrarily weight the angular distribution of the primary lepton to obtain an approximate 
agreement with the data. Although, as shown in Appendix B, the remaining kinematics is 
poorly modeled by the hypothesis of a final state containing two massive particles ( W and 
N), this simulation provides a jet multiplicity distribution consistent with the simpler N 
production. A fit of the rates of tagged W+ jet events with the SM simulation implemented 
with this WN production also returns a number of tagged WN events consistent with the 
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fit shown in Tables [TJ and |T|. The detector acceptance for such a WN process is about 9% 
without including the W leptonic branching ratio. 

In a second attempt to evaluate the detector acceptance, we model the production of 
events with a superjet with a more general 2 — > 4 hard scattering process. As detailed 
in Appendix C, the matrix element of this process is tuned using an effective Lagrangian 
approach to account for the salient kinematic features of the data. In this picture the 
primary lepton and the missing energy do not come from W decays. In this simulation, the 
yield of events before and after tagging as a function of the jet multiplicity is similar to that 
of the pp N simulation (Tables | to |TTTD . Under the assumption that the lepton in the 



final state can be an electron, a muon or a r with equal probability, the detector acceptance 
for this simulated process is 10.8%, consistent with the WN simulation. 

If we average the two techniques for estimating the detector acceptance, then the 52.8 
± 22.1 events with a superjet required by the fit in Table |I| correspond to a cross section of 
about 5 ± 2 pb. 

IV. CONCLUSIONS 

We present additional studies of the kinematics of the anomalous events with a superjet 
observed by CDF. We model several properties of the superjets with the hypothesis of the 
production and decay of a light scalar quark. Using phenomenological approaches, based 
on this ansatz, to evaluate the detector acceptance and the tagging efficiencies for events 
with a superjet, we estimate that their production cross section is approximately 5 ± 2 pb, 
independent of the details of the kinematic distributions. At the same time, we derive a new 
estimate of the tt cross section value, 4.0 ±1.5 pb, using W+ > 1 jet events with SECVTX, 



JPB and SLT tags. This result agrees with the D0 measurement (TT|, o t i = 4.1 ± 2.0 
pb, derived in the same topological channel by using kinematical cuts designed to select tt 
events. 
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APPENDIX A: SIMULATION OF THE PROCESS pp -» N 



We simulate the process pp —> N with the herwig Monte Carlo program |12| . We 
produce the state N using option 1605, which calculates the production of a Higgs boson in 
pp interactions. The Higgs boson is forced to decay to bb. The decay of the 6-hadron formed 



by herwig is performed using the QQ Monte Carlo program [13]. The 6-quark is transformed 
into a scalar quark in the following way. Having modified the herwig generator, we set the 
6-quark mass to 3.6 GeV/c 2 . The 6-quark is hadronized by herwig as a fermion. At the 
end of the hadronization process, when a 6-hadron has been formed, we change its identity 
to a fictitious 6-hadron, but we keep the hadron mass evaluated by herwig at this stage. 
The decay of the 6-hadron is modeled with herwig using the spectator model (routine 
hwdhvy). In this model, herwig weights in the routine hwdhqk the phase space of the 
three-body semileptonic decay of a 6-quark with the V-A matrix element calculated with 
the routine hwdwwt. In the case of a scalar quark we modified herwig to weight the 
phase-space with the matrix element: 

— — = K[(l - z c )(l - zi) -R, + R c (z c - Zl + Ri>- R c )} (Al) 

dz c dzi 

where K is a normalization factor, R c = m 2 c /m^ and R„ = m?/m?. 
Here z c and z t are defined as 

2Pb • Pc , 2p~ b ■ pi 
z c = 5— and zi = 5— 

mf mi 

b b 
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The phase space limits are: 



1\ R C < z c < 1 + R c - Ry 



1 + R c — R v — z c 1 + R c — R u — z c 
< zi < 



4i? c ]/2 



\Z C + \ Z, 



4R c }/2 



This matrix element for the decay b — > civ is derived from the tree level calculation outlined 



in Ref. 14 



M 



y -a 2 v n F L \v ci \ 

j=l J 



Ujicosd ■ 



m b Uj2sin9 
\/2myv cos (3 



m c V*2 cos Q 

m Xj u{pc)P R v(pi) + — —u(Pc)^Prv(j>i) 

V 2m\Y sin p 



(A2) 



where x/ s are the chargino mass eigenstates, and U and V* are the mixing matrices for 
the right and left-handed charginos, respectively. The first subscript of U, V corresponds 
to mass eigenstates and the second to weak eigenstates (1 for the gaugino and 2 for the 
higgsino). Here tan/? is the ratio of the vacuum expectation values of the two Higgs fields, 9 
is the mixing angle between left-handed and right-handed scalar quarks, is the fraction 
of left-handed component of the scalar neutrino, and \V c i\ is the CKM matrix element. 
Equation (Al) follows from (A2) when the decay is mediated by the higgsino coupling to 
the right-handed matter. If the decay is mediated by the gaugino coupling to the left-handed 
matter the matrix element is 



dY 



K(z c + Zl -l + R v -R c ). 



(A3) 



dz c dzi 

In the latter case the two fermions in the final state are both left-handed and tend to 
be produced back-to-back since the initial state is spinless. Using equation (Al) the two 
fermions in the final state have opposite handedness and tend to be produced more collinear 
than when using equation (A3). It follows that the matrix element (Al) produces leptons 
with a momentum distribution appreciably softer than when using the matrix element (A3) 
or a phase space decay, as done in Ref. || . The c-quark emerging from the 6-quark decay is 
then recombined in herwig with the spectator quark by the routine hwcfor. The decay 
of the resulting c-hadron is performed with the QQ Monte Carlo program. In the spectator 
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model, excited D-meson states are produced only by the hadronic current carried by the 
virtual W (so called upper vertex). Since we impose that the gauge fermion involved in the 
scalar quark decay has only leptonic decay modes, the simulation produces a very simple 
list of c-hadrons with respect to the QQ generator (see Figure [|). 



10' 



10 



p 

3 

a io 2 
■— 



10 



10' 



io- 



p 

3 

5 10 2 

M 

■4—1 



1.5 



5.5 



b-hadrons 



6.5 



M (GeV/c ) 



c-hadrons 
from b decay ^ 



II 



10' 



10 1 r 



P 
3 

a io 2 : 



10 



as 



10' 



I0 1 r 



P 
3 

a io 2 : 

■I— I 



2.5 



3 



3.5 



1.5 



3.5 



b-hadrons 



4.5 



M (GeV/c ) 



c-hadrons 
from b decay _ 



2.5 



3.5 



M (GeV/c ) M (GeV/c ) 

FIG. 3. Mass spectrum of the b and 6-hadrons produced by the herwig simulation using 
rrib=A.75 GeV/c 2 (a) and m^=3.65 GeV/c 2 (b). (c): mass spectrum of the c-hadrons from 6-hadron 
decays modeled with QQ. (d): mass spectrum of the c-hadrons from 6-hadron decays modeled with 
the herwig spectator model. 
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APPENDIX B: SIMULATION OF THE PROCESS pp -» WN 

We generate the process q^qj — > WH with m H = 220 GeV/c 2 using option 2605 of 



herwig and the MRS(G) set of structure functions fig] . We identify the Higgs boson with 
the state N and we force it to decay to 66. As shown in Figure |J the distributions of 
the pseudo-rapidity and transverse momentum of primary leptons in this simulation and 
in the data are quite different. We attempt to model the data by considering the process 
qiqj — ► WN as a 2 — ► 3 hard scattering (g^- — > Z + El T iV, with subsequent decay iV — > 66). 
For the sake of simplicity, we ignore the polarization of the outcoming W calculated by 
herwig and decay the W boson into Iv according to two-body phase space. Then, in order 
to model the observed pseudo-rapidity distribution of the primary leptons, we weight their 
distribution in the q^j rest-frame with the function z A = cos 4 ^, where Q\ is the polar angle 
of the primary lepton with respect to the gj-dfrection. As shown in Figure H this attempt 
is quite successful. 

Using this simulation, we compare additional kinematical quantities to the data. Figure |6| 
shows the correlation between the pseudo-rapidity of the primary lepton and the rapidity 
of the di-jet system N. In the simulation, in which the W boson recoils against the state 
N, primary leptons tend to have rapidities opposite to the N direction. In the data, the N 
rapidity is quite central, independent of the primary lepton pseudo-rapidity. In the data, 
primary leptons are produced at large rapidities while the state N is produced quite centrally; 
therefore, most of the N transverse momentum is balanced by the missing transverse energy 
(Figure 0). 

Pseudo-rapidity distributions of the 6-jets in the data and in the WN simulation are 
shown in Figure || Differently from the small sample of data, simulated 6-jets produced 
by the decay of a massive state N tend to fill quite uniformly the 7/-region covered by the 
detector. 
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FIG. 4. Transverse momentum (a) and pseudo-rapidity (b) distributions of the primary leptons 

resulting from the decay of W bosons produced in association with a Higgs of mass 220 GeV/c 2 ; 
(c) and (d) are the analogous distributions for the data. 
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FIG. 5. Transverse momentum and pseudo-rapidity distributions of the primary lepton. The 
data (•) are compared to a q^qj — > /$ T iV simulation in which the angular distribution of the primary 
lepton is weighted with the function cos 4 9i (see text). 
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data (a) and in the modified WN simulation (b). The distribution of the pseudo-rapidity f] su j of 
the superjets versus r\\ is shown in (c) for the data and in (d) for the simulation. 
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FIG. 7. Distribution of # T vs p!f in the 13 events with a superjet. 
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APPENDIX C: EFFECTIVE LAGRANGIAN APPROACH 

The matrix element used in Sec. Ill has been derived using an effective Lagrangian 
approach to incorporate the salient features of the data. The production of events with a 
superjet is modeled with the 2 — > 4 hard scattering process 

u( Pl ) + d(p 2 ) -> e + (pi) + u s {k) + b( Pb ) + b s (p s ) (CI) 

where b s and u s are a scalar quark and a scalar neutrino. A scalar neutrino is used in the 
final state to model the measured missing energy for no compelling reason other than the 
fact that we have already identified superjets with a scalar quark. 
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One mechanism capable of producing a angular distribution for the primary leptons 
in the final state is a high derivative coupling to the lepton and quark fields. As a first step, 
using the scalar fields <fi of the v s and (f) of the b s , the following dimension-1 and spin-2 fields 
are constructed: 

</v(x) = (1/A)%d u <j)(x) and 0^ = {1/A)%d v ${x) 

where A sets the mass scale of the interaction. In analogy, from the lepton and quark spinors 
one constructs the following dimension-3/2 and spin-5/2 fields: 

(x) = {l/Afd^ipiix); O) = {l/Afd^d v ip h {x) 

and 

= (1/A) + ±^l> d (x)] 

= (1/A) W[±^»C«0 + l^M*)] 

where s q is the sign of the initial state (u + d) charge. The peculiar forms Xq V an d f i 
used to account for the fields of the initial state partons, are an attempt to also model the 
asymmetries observed in the rapidity distributions of the data. 

These fields are used to write the following high- derivative, but point-like, effective La- 
grangian: 

= (^)dM^)^( x ^T(^)}Yd x {^ WW d u [(drC s Mhp K$» ( d °xf (*))]} 

where A du B = A(d u B) — (d u A)B, and ip" is the transpose of ip' , the "charge-conjugate" 
of ijj. To avoid conflict with unitarity in the high energy limit, we introduce a form factor 
to compensate for the twenty derivatives used above []. However, the effective Lagrangian 

1 Since we use derivative couplings to fundamental fields, it is natural to assume that the hypo- 
thetical object involved in the s-channel exchange is composite. 
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acquires a reasonable s-behaviour only after the inclusion of an additional s-channel prop- 
agator. The squared amplitude for the process (CI) (and its charge conjugate), averaged 
over initial quark colors, then reads 



2/-:, .... AE, 



1 ^ 2 1- ^ 



. l-cos^ fc l + cos# a 
coswz smWb ( )( 



In the kinematics of Eq. (C2), all six particles in the initial and final states are taken to be 
massless. The energies and polar angles [ (Ei,$i) for the lepton, (E b ,-d b ) for the b-quark, 
(k,$k) for the scalar neutrino, and (E s ,i9 s ) for the scalar b-quark] are defined in the ud 
center of mass system, and u(pi) is taken to be the in the z-direction. Such an effective 
Lagrangian has been constructed by iterations in which we compared data and simulation. 
We have started using a dimension-6 Lagrangian and we kept increasing the Lagrangian 
dimension till we obtained an acceptable description of the data. 



The process (CI) has been kindly implemented into the herwig generator |IB |. In 
the calculation of the matrix element (C2), we use / = 1 and Ai = GeV. In order to 
approximately reproduce the distribution of the invariant mass of the final state in the data, 
we use M ~ 350 GeV/c 2 and T ~ 5 GeV/c 2 |. The cross section value of 5 pb derived in 
Sec. Ill corresponds to an interaction scale A of approximately 10 MeV. 

Figures |9] to [II] compare pseudo-rapidity, transverse momentum and invariant mass dis- 
tributions in the data and the simulation normalized to the same number of events. Data 
and simulation are in qualitatively good agreement. 



2 Widths r as large as 50 GeV/c 2 and masses M in the range of 300 — 400 GeV/c 2 also provide 
quite similar results. 
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FIG. 9. Transverse momentum and pseudo-rapidity distributions of primary leptons and super- 
jets. 
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FIG. 10. Transverse momentum and rapidity distributions of 6-jets and of the system N 
sisting of the 6-jet and the superjet. 
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FIG. 11. Distributions of: Ei T (a); the transverse momentum of the system l$ T composed by 
the primary lepton and the transverse missing energy (b); the transverse mass of the system l$ T 
(c); and the invariant mass of the system ./V consisting of the 6-jet and the superjet (d). 
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